Elementary events of electron transfer in a voltage-driven quantum point contact 
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We find that the statistics of electron transfer in a coherent quantum point contact driven by 
an arbitrary time-dependent voltage is composed of elementary events of two kinds: unidirectional 
one-electron transfers determining the average current and bidirectional two-electron processes con- 
tributing to the noise only. This result pertains at vanishing temperature while the extended 
Keldysh-Green's function formalism in use also enables the systematic calculation of the higher- 
order current correlators at finite temperatures. 
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The most detailed description of the charge transfer in 
coherent conductors is a statistical one. At constant bias, 
the full counting statistics (FCS) of electron transfer [lj 
can be directly interpreted in terms of elementary events 
independent at different energies. The FCS approach 
is readily generalized to the case of a time-dependent 
voltage bias @, 0] . The current fluctuations in coherent 
systems driven by a periodic voltage strongly depend on 
the shape of the driving [3] , which frequently is not ap- 
parent in the average current The noise power, for 
instance, exhibits at low temperatures a piecewise linear 
dependence on the dc voltage with kinks corresponding 
to integer multiples of the ac drive frequency and slopes 
which depend on the shape and the amplitude of the ac 
component. This dependence has been observed experi- 
mentally in normal coherent conductors 6 1 and diffusive 
normal metal-superconductor junctions 7|. 

The elementary events of charge transfer driven by a 
general time-dependent voltage have not been identified 
so far. The time dependence mixes the electron states at 
different energies [5| which makes this question both in- 
teresting and non-trivial. The first step in this research 
has been made in Q for a special choice of the time- 
dependent voltage. The authors have considered a su- 
perposition of overlapping Lorentzian pulses of the same 
sign ("solitons"), with each pulse carrying a single charge 
quantum. The resulting charge transfer is unidirectional 
with a binomial distribution of transmitted charges. The 
number of attempts per unit time for quasiparticles to 
traverse the junction is given by the dc component of the 
voltage, independent of the overlap between the pulses 
and their duration Q. It has been shown that such su- 
perposition minimizes the noise reducing it to that of a 
corresponding dc bias. A microscopic picture behind the 
soliton pulses has been revealed only recently (lfj |. In 
contrast to a general voltage pulse which can in princi- 
ple create a random number of electron-hole pairs with 
random directions, a soliton pulse at zero temperature 
always creates a single electron-hole pair with quasipar- 
ticles moving in opposite directions. One of the quasipar- 



ticles (say, electron) comes to the contact and takes part 
in the transport while the hole goes away. Therefore, 
soliton pulses can be used to create minimal excitation 
states with " pure" electrons excited from the filled Fermi 
sea and no holes left below. The existence of such states 
can be probed by noise measurements 1^, U , l2| • 

In this Letter, we do identify the independent elemen- 
tary events for an arbitrary time-dependent driving ap- 
plied to a generic conductor. Since generic conductor at 
low energies can be represented as a collection of indepen- 
dent transport channels, it is enough to specify elemen- 
tary events for a single channel of transmission T. The 
answer is surprisingly simple. There are two kinds of such 
events: We call them bidirectional and unidirectional. In 
the course of a bidirectional event k an electron-hole pair 
is created with probability sin 2 (afc/2), with ctk being de- 
termined by the details of the time-dependent voltage. 
The electron and hole move in the same direction reach- 
ing the scatterer. The charge transfer occurs if the elec- 
tron is transmitted and the hole is reflected, or vice versa 
[Fig. QJa,b)]. The probabilities of both outcomes, TR 
(R being reflection coefficient), are the same. Therefore, 
the bidirectional events do not contribute to the average 
current and odd cumulants of the charge transferred al- 
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FIG. 1: Schematic representation of elementary events: bidi- 
rectional (a, b) and unidirectional (c). Shifts of the effective 
chemical potential in the left lead due to time-dependent volt- 
age drive are indicated by shading. For periodic drive, the dc 
voltage component [panel (d) , dash-dotted line] describes uni- 
directional charge transfer, while the ac component (dashed 
curve) describes bidirectional events affecting the noise and 
higher-order even cumulants. 
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though they do contribute to the noise and higher-order 
even cumulants. A specific example of a bidirectional 
event for a soliton-antisoliton pulse was given in Q . 

The unidirectional events are the same as for a con- 
stant bias or a soliton pulse. They are characterized by 
chirality Ki = ±1 which gives the direction of the charge 
transfer. An electron-hole pair is always created in the 
course of the event, with electron and hole moving in op- 
posite directions [Fig. QJc)] . Either electron (k; = 1) or 
hole (ki = —1) passes the contact with probability T, 
thus contributing to the current. 

Mathematically, the above description corresponds to 
the cumulant generating function S(x) — Si(x) + S 2 (x)i 
where 
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are the contributions of the bidirectional and unidirec- 
tional events, respectively. Here x is the counting field, 
and afc and k\ are the parameters of the driving to be 
specified later. The sum in both formulas is over the set 
of corresponding events [l3T ] . The elementary events have 
been inferred from the form of the cumulant-generating 
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function, as it has been done in ^ 

The cumulant-generating function given by Eqs. ([T]) 
and ([2]), together with the interpretation, is the main re- 
sult of this Letter. It holds at zero temperature only: 
since the elementary events are the electron-hole pairs 
created by the applied voltage, the presence of thermally 
excited pairs will smear the picture. Equations (fTJ) and 
(|2|) contain the complete %-field dependence in explicit 
form which allows for the calculation of higher-order cu- 
mulants and charge transfer statistics for arbitrary time- 
dependent voltage. The probability that N charges are 
transmitted within the time of measurement is given by 
P(N) = (2TT)- 1 f^dxexp[S(x) - iNx]- Higher-order 
derivatives of S with respect to x are proportional to 
the cumulants of transmitted charge, or equivalently, to 
higher-order current correlators at zero frequency. The 
details of the driving are contained in the set of param- 
eters {afc} and separated from the x-field dependence. 
This opens an interesting possibility to excite the spe- 
cific elementary processes and design the charge transfer 
statistics by appropriate time dependence of the applied 
voltage, with possible applications in production and de- 
tection of the many-body entangled states [1 

Below we present the microscopic derivation of Eqs. 
and |2]). We neglect charging effects and assume in- 
stantaneous scattering at the contact with quasiparti- 
cle dwell times much smaller than the characteristic 
time scale of the voltage variations. The approach we 



use is the nonequilibrium Keldysh-Green's function tech- 
nique, extended to access the full counting statistics 
m El H ED]. The Green's functions of the left (1) 
and right (2) leads are given by 2(| 21 1 



G x = e-^/ 2 (I 2 _ty e}^l\ G 2 = (l , (3) 

where fi = ( 5 q ) is a matrix in Keldysh(") space. Here- 
after we use a compact operator notation in which the 
time (or energy) indices are suppressed and the products 
are interpreted in terms of convolution over internal in- 
dices, e.g., {G 1 G 2 )(t',t") = J dhG^t' M)G 2 {h,t") (and 
similar in the energy representation). The equilibrium 
Green's function G-2.it' — t") depends only on time differ- 
ence. In the energy representation 62(8', £") is diagonal 
in energy indices with h(£ £") = tanh(£'/2T e ) 2ir5{£' - 
£"). Here the quasiparticle energy £ is measured with 
respect to the chemical potential in the absence of the 
bias and T e is the temperature. The Green's func- 
tion Gi(t',t") depends on two time (or energy) argu- 
ments. It takes into account the effect of applied volt- 
age V(t) across the junction through the gauge trans- 
formation h = UhU^ which makes G\ nondiagonal in 
energy representation. The unitary operator U is given 
by U(t',t") = f(t')5(t' - t") in the time representation, 

where f(t') = exp[— i J eV(t)dt\. The cumulant gener- 
ating function S(x) of the charge transfer through the 
junction is given by [2l], EH 



S(x) = Trln 
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Here the trace and products of Green's functions include 
both summation in Keldysh indices and integration over 
time (energy). For a dc voltage bias, G\ and G 2 are 
diagonal in energy indices and S(x) is readily interpreted 
in terms of elementary events independent at different 
energies 21] . To deduce the elementary events in the 
presence of time dependent voltage drive it is necessary to 
diagonalize {Gri, G 2 } £l£ „. The diagonalization procedure 
is described in the following. 

For the anticommutator of the Green's functions we 
find {G 1 ,G 2 }/2 - 1 = -2sin(x/2)(.A + B), with A = 
(Jg) (g) A, and B = (g^ 6 ) ® B. Here A = (1 - 

Wi)sin(x/2) + i(/i-ft,)cos(x/2), B = (1 — hh) sin(x/2) + 
i(h — /i)cos(x/2), b = — icot(x/2), and £g> is the ten- 
sor product. Since AB — BA — 0, the operators A 
and B commute and satisfy for integer n: (A + B) n = 
( h 0) ® An + (0 ~i) ® Bn - Therefore, S(x) given by 
Eq. ([2]) reduces to 



S = Trln 
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A further simplification of S(x) is possible in the zero 
temperature limit, in which the hermitian /i-operators 
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are involutive: h 2 = h 2 = 1. The operators hh and hh 
are mutually inverse and commute with each other. Be- 
cause hh is unitary, it has the eigenvalues of the form 
e lak with real a k , and possesses an orthonormal eigenba- 
sis {v Qfc }. The typical eigenvalues of hh (or hh) appear 
in pairs e ±M with the corresponding eigenvectors v a and 
v_ a = hv a . In the basis (v Q ,v_ Q ) operators hh and 
hh are diagonal and given by hh — diag(e M ,e~ M ) and 
hh = diag(e~* a , e la ). The eigensubspaces span(v a ,v_ Q ) 
of the anticommutator {h, h} are invariant with respect 
to h, h, and A because of [h, {h, h}] = [h, {h, h}] = 0. 
The operators h and h are anti-diagonal in the basis 
(v Q ,v_ Q ), with matrix components hyi = 1, h 2 i = 1, 
hi2 = e~ la , and h 2 i — e la . The operator A can be diag- 
onalized in invariant subspaces, with typical eigenvalues 
given by 

ev A = 2 sin(a/2) ( sin(a/2) sin(x/2) 

±i^l-sin 2 (a/2)sin 2 (x/2)). (6) 

Similarly, we obtain ev B = ev A. From Eqs. (J3J) and © 
we recover the generating function 5*1 given by Eq. |T]), 
which is associated with the paired eigenvalues e ±lQfc . 

There are, however, some special eigenvectors of hh 
which do not appear in pairs. The pair property dis- 
cussed above was based on the assumption that v a and 
hv a — v_ Q are linearly independent vectors. In the spe- 
cial case, these vectors are the same apart from a co- 
efficient. Therefore, the special eigenvectors of hh are 
the eigenvectors of both h and h with eigenvalues ±1. 
This means that the special eigenvectors posses chirality, 
with positive (negative) chirality defined by hv — v and 
hv = —v (hv = —v and hv = v). From Eq. (|5|) we obtain 
the generating function ^(x) given by Eq. ([2]), where I 
labels the special eigenvectors and ki is the chirality. 

In the following we focus on a periodic driving V(t + 
t) = V(t) with the period r = 2tt/lo, for which the 
eigenvalues of hh can be easily obtained by matrix diag- 
onalization. The operator h couples only energies which 
differ by an integer multiple of lo, which allows to map 
the problem into the energy interval < £ < lo while 
retaining the discrete matrix structure in steps of lo. 
Therefore, the trace operation in Eq. ^ becomes an 
integral over £ and the trace in discrete matrix indices. 
The operator hh in the energy representation is given by 
(hh) nm {£) = sign(£ + nw) fn+kf^+k si S n ( £ ~ kuj ~ 
eV), with f n = (1/t) f_^ /2 dt e- l ti dt ' ^^'''e"'. Here 
V = (1/r) J V(t)dt is the dc voltage offset and AV(t) = 
V(t) — V is the ac voltage component. The coefficients 
/„ satisfy J^k fn+kfm+k = 5 «™ and T, n n \fn\ 2 = 0. 

To evaluate S(x) for a given periodic voltage drive V(t) 
it is necessary to diagonalize (hh) nm (£). First we ana- 
lyze the contribution of typical eigenvalues e ±m . The 
matrix (hh) nm {£) is piecewise constant for £ G (0,wi) 



and £ G (lj%,u)), where lj% = eV — Ncj and N — [eV/u\ 
is the largest integer less than or equal eV /lo. The eigen- 
values e ±lQfci ( R > of (hh) nm are calculated for £ 6 (0,u>i) 
[£ G (u)i,u))] using finite-dimensional matrices, with the 
cutoff in indices n and m being much larger than the char- 
acteristic scale on which |/„| vanish. Further increase 
of the size of matrix just brings more eigenvalues with 
ctfc = which do not contribute to 5*(x), and does not 
change the rest with a k 0. This is a signature that 
all important Fourier components of the drive are taken 
into account. The eigenvalues e ±mfcL < R > give rise to two 
terms, Si = Sil + Sir, with 

SilMx) = M L , R Y / Hl+TRsm 2 (a k L,R/2) 

k 

x (e ix + e~ ix - 2)}. (7) 

Here Ml = toLUi/ir, Mr = to(u> — u>i)/tt, and to is the 
total measurement time which is much larger than t and 
the characteristic time scale on which the current fluctu- 
ations are correlated. 

The special eigenvectors all have the same chirality 
which is given by the sign of the dc offset V. For 
eV > 0, there are Ni = N + 1 special eigenvectors 
for £ £ (0, u>i) and N 2 = N for £ G (uii,u). Because 
eV — Niwi + N 2 (uj — uji), the effect of the special eigen- 
vectors is the same as of the dc bias 

S 2 (X) = — ln[l+T(e- <x -l)]. (8) 

TT 

Comparing Eqs. and (jSJ) we see that unidirectional 
events for periodic drive are uncountable. The summa- 
tion in Eq. ([2]) stands both for the energy integration in 
the interval lo and the trace in the discrete matrix in- 
dices. In the limit of a single pulse lo — > unidirectional 
events remain uncountable for a generic voltage, while 
being countable, e.g., for soliton pulses carrying integer 
number of charge quanta Q . 

Equations (J7J) and ||5J) determine the charge trans- 
fer statistics at zero temperature for an arbitrary peri- 
odic voltage applied. The generating function consists 
of a binomial part (S2) which depends on the dc off- 
set V only, and a contribution of the ac voltage com- 
ponent (Si) [Fig. [ljd)]. The latter is the sum of two 
terms which depend on the number of unidirectional at- 
tempts per period eV /lo. The simplest statistics is ob- 
tained for an integer number of attempts for which Sil 
vanishes 0]. The Fourier components of the optimal 
Lorentzian pulses Vj,(f) = (2T L /e) ^2 k [{t — kr) 2 + t 2 ]^ 1 
of width tl > are given by /_i = — e - 27rr i7' r j f n = 

e -2nnr L /r _ e -2n(n+2)r L / r for n > Q, and /„ = Oth- 
erwise. In this case Sil = Sir = and the statistics 
is exactly binomial with one electron-hole excitation per 
period, in agreement with Refs. (jl. Hoi]. 

The elementary events at zero temperature can be 
probed by noise measurements. For example, in the case 
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charge transfer is the simplest for an integer number of 
attempts for quasiparticles to traverse the junction. This 
results in photon-assisted effects in even-order cumulants 
as a function of a dc offset. The approach we have used 
also allows for the systematic calculation of higher-order 
cumulants at finite temperatures. 
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tov and C. Bruder. This work has been supported by 
the Swiss NSF and NCCR " Nanoscience" (MV), the 
DFG through SFB 513 and the Landesstiftung Baden- 
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FIG. 2: The probability of elementary events for harmonic 
drive with amplitude Vo (upper panel). With increasing am- 
plitude more and more eigenvalues come into play and 
contribute to transport. The derivative of the noise power 

with respect to Vo ( — ) decomposed into contributions ( ) 

of elementary events (lower panel) . 

of an ac drive with V = 0, only bidirectional events of 
i?-type remain [S(x) — Sm(x)]- Both the number of 
events and their probabilities increase with increasing the 
driving amplitude Vq, which results in the characteristic 
oscillatory change of the slope of the current noise power 
Pi = (4e 2 w/7r)T(l - T) J2 k shr 2 (a fc /2). The decomposi- 
tion of dPi/dVa into contributions of elementary events 
for harmonic drive is shown in Fig. O 

Our method also enables the efficient and systematic 
analytic calculation of the higher-order cumulants at fi- 
nite temperatures. They can be obtained directly from 
Eq. by expansion in the counting field to the certain 
order before taking the trace. The trace of a finite num- 
ber of terms can be taken in the original basis in which 
Gi and G*2 are defined. The details of this approach will 
be given elsewhere. However, the formulas obtained (as 
a function of {f n }) can not be interpreted as elementary 
events term by term. To identify the elementary events 
it is necessary to find S(x) which requires full expansion 
or diagonalization, as presented above. 

In conclusion, we have studied the statistics of the 
charge transfer in a quantum point contact driven by 
time-dependent voltage. We have deduced the elemen- 
tary transport processes at zero temperature from an 
analytical result for the cumulant generating function. 
The transport consists of unidirectional and bidirectional 
charge transfer events which can be interpreted in terms 
of electrons and holes which move in opposite and the 
same directions, respectively. Unidirectional events ac- 
count for the net charge transfer and are described by 
binomial cumulant generating function which depends on 
the dc voltage offset. Bidirectional events contribute only 
to even cumulants of charge transfer at zero temperature. 
They are created with probability which depends on the 
shape of the ac voltage component. The statistics of 
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